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SUMMARY 

Steady, diabatic (nonadlabatic) , frictional, vafiable-area flOT 
of a compressible fluid is treated in differential form on tbe .basis 
of the one -dimensional approximation; Tbe basic eguations are first 
stated in te 2 rms of pressure, temperature, density, and velocity of 
the fluid. Considerable simplification and unification of the eq.ua- 
tions is then achieTed by choosing tbe square of the local Mach num- 
ber as one of the variables to describe the flow. 

The transformed system of equations thus obtained is first 
examined with regard to the existence of a solution. It is shown 
that, in general, a solution exists -vdiose calculation requires know- 
ledge only of 'the variation with position of any three of the 
dependent variables of tbe system. The direction of change of the 
flow variables oan be obtained directly from tbe transformed eqiia- 
tlons without integration. As examples of this application of the 
equations, the direction of change of the flow variables is deter- 
mined for two special flows. 

In the particular case when the local toch number M = 1, a 
special condition must be satisfied by the flow if a solution is to 
exist. This condition restricts tbe Joint rate of variation of 
heating, friction, and area at M = 1. Further analysis indicates 
that when a solution exists at this point it is not necessarily 
unique . 

Finally it is shown that the physical phenomenon of choking, 
which is known to occur in certain simple flow situations, is 
related to restrictions imposed on the variables by the form^of the 
transformed equations. The phenomenon of choking is thus given a 
more general significance in that tbe transformed equations apply 
to a more general type of flow than has hitherto been treated. 
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INTRODUCTION 

The rational experimental development of Jet- eind rocket- 
propulsion power plants requires adequate knowledge of tlje theoreti- 
cal mechanics of diabatic (nonadiabatic) , frictional, variable-area 
compressible fluid flow. The differential equations describing this 
typo of flow are well known. (See, for oxamplo, references 1(a), 
1(b), 2, 3, and 4.) Their solution in the three-dimensional case, 
however, is so difficult that some simplification is necessary to 
permit development of the theory in a form immediately useful for 
technical applications. 

In the present paper, such simpliflcatlan is effected by gen- 
eralizing the familiar "one -dimensional" or hydraulic treatment of 
fluid flow to include the simultaneous effects of heat addition, 
friction, and area change upon the flow of a compresaiblo fluid 
rather than by attempting to show that the one -dimensinml -approx- 
imation follows from a simplification of the hydrodynamic and heat- 
flow eq'uations in their general three-dimensional form. The gen- 
eralization leads to one -dimensional equations in differential form, 
which are identical with equations previously used by other inves- 
tigators in less general cases. 

Generalized conservation equations have been derived in appeai- 
dix A in order that a complete and logical basis for the theory may 
be accessible to the reader. The resulting theory is intended to 
serve as a foundation in differential form for calculatxon of ail 
types of mathematically continuous (that is, shocklcss) flow of per- 
fect gases to which the one-dimensional approximation is applicable. 
Thus the theory applies directly to compressible flow in combustion 
chambers and also, with but slight modification, to flow in turbines 
eind compressors (cf. reference 5) and nozzles and diffusers whenever 
the one- dimensional approximation is valid. 

In order to obtain convenient and unified equations, the gen- 
eralized relations are transformed by introducing a new easlc 
variable, the square of the local 11a ch number ^ N. Pressure 
and temperature are chosen as the additional basic variables; other 
relevant flow variables (for example, density, velocity, mass flow) 
may be expressed in terms of Mach number squared, pressure, and tem- 
perature. Values of M from zero to Infinity are considered; the 
treatment is therefore applicable to- both subsonic and supersonic 
flow. 



NACA TN Wo, 1336 


3 


The variable M has been used throughout dlTferential treat- 
ments by Gukhman (reference 6), Bailey (reference 7); Nielsen 
(reference 8), who investigated various examples of frictional, 
diabatic compressible flow. A related variable 

Z = (7-1)m2/[8 + (7-1)m 2], which can be used alternatively with M, 
will be discussed briefly in appendix B. Pertinent papers in which 
M is not used extensively are references 9 and 10, which report studies 
of Isothermal etnd of adiabatic frictional flow, respectively. A 
treatment of frictionless diabatic compressible flow carried out by 
Szczeniowski (reference-11) is partly in differential form. The 
same sub;1ect, using the M language without differential formula- 
tion, is discussed in reference 12. The variable M has also bean 
employed to advantage in reference 5 for analysis of cwnpressible 
flow through turbines and compressors, a related field that is not 
specifically discussed in the present paper. 

In the general case, the differential eqiiations obtained in the 
present treatment do not permit of formal integration: but being of ^ 
the first order, they axe particularly amenable to numerical methods. 

A solution of the system is shown to exist, except possibly at sonic 
velocity, and the behavior of the solution in this neighborhood is 
investigated. From the differential eq.uatlons viseful. information 
may be easily obtained about direction of changes in the flow vari- 
ables, Choking is shown to be a consequence of a certain property 
of the equations. 


THE ONE-DIMEINSIONAL APPROXIMATION 
Basic Equations 

The "one -dimensional" steady-flow theory utilizes a model con- 
sisting of a perfect gas contained within a duct, across any section 
of which the flow variables are constant. Only the component of 
velocity normal to the section is considered; body forces are neg- 
lected, and heat, whether supplied by combustion, conversion of 
frictional work, or conduction from the walls, is assumed to be 
transferred insteintaneously and completely bixt only transversely 
throughout the cross section, which may be of variable area. Each 
flow variable can thus be considered as a function of a single 
parameter, say the distance along the axis of the tube, whence the 
term "one dimensional," 
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The conventional variables — pressure, temperature, density, 
and velocity in one -dimensional flow — are connected by four rola- 
tions derivable from the first law of thermodynamics, the conserva- 
tion of mass, the second law of motion, and the thermal oauation of 
state for a perfect gas. 

The four relations are: 


Conservation of energy 

Cp dT + VdY = dQ 

(1) 

Conservation of mass 

d(pVA) = 0 

(2) 

Eqiiation of motion 

~dp a pVdV + pdF 

(3) 

Equation of state 

d(p/ElpT) = 0 

(4) 


The specific heat- at constant pressure Cp and the gas constant E 
do not vary in the flow. The symbols p, T, p, and T, respec- 
tively, stand for density, velocity, absolute statl-c pressure, and 
absolute static temperature. The pipe area, which may bo vazdablo, 
is represented by A. Heat added per \mit mass is Indicated by Q, 
and work per unit mass done against friction by E. Consistent 
units are used throughout-. In ecLuations (1) tbroiagh (4) each vari- 
able* is to be considered as a functicsn of a single paxeunetor, such 
as the distance x along the tube considered positive in the direc- 
tion of flowj and, of course, the meaning of each differential du 
is then given by 


du = u' (x)dx 

Equations (1) to (3) are customarily used without explicit 
recoij^ition of their true meaning with regard to the one -dimensional 
approximation. The interpretation of the quantity dF in particular 
is often obscTxre. In order to provide a logical, unified basis for 
the theory, equations (l) to (3) are derived in appendix A; spoclal 
care is taken to keep the dorlvations within the framework of the 
one-dimenalonal approximation. “t°- 


Applicability 

The validity of the onc-dimensional apuroxlmation depends uijon 
the assumption of the ;iniformity of flow condltdons across a nlano 
normal to the direction of flow. Experience has shown that this 
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assumption constitutes an adequate approximation in many special 
cases; in particular, with subsonic turbulent flow in long pipes 
without separation, the .reasonably flat velocity profile permits 
the use of equations derived on this basis. In cases involving 
incomplete growth of boundary layer or where separation of flow occurs, 
however, there is grave doubt as to the applicability of a one- 
dimensional treatment. Although bomidary- layer effects are somewhat 
amenable to calculation, the occurrence of separation is difficult 
or impossible to predict and the question of applicability must usually 
be determined by experiment or estimated by experience. 

The one -dimensional approximation would not be valid if oblique 
shocks occur in the flow. Nor can normal shocks, if treated as flow 
discontinuities, be bandied in the differential form of the present 
approximation. If, however, in equations (l) and (3), dQ and dP 
are considered to depend upon the derivatives of T and T and if 
heat and momentum transfer in the direction of flow is allowed, then 
the equations for continuous norasal shock (reference 1(c), p. E19) 
can be put in the form of equations (1) to (4). 

In the development and use of equations (1) to (4) various 
approximations are made, such as neglecting the squares of velocity 
components normal to the direction of flow, replacing the square of 
the cosine of the half -angle by unity, and assuming the constancy 
of E and Cp. In this paper no attempt is made to state under 
what circumstances such approximations are suitable. 


TRANSFOBMATION OF EQUATIONS 
Cheinge of Variables 

A canonical form for equations (1) to (4) is obtained by taking 
logarithmic derivatives and choosing as a variable the square of the 
local Mach number 


M^= N = Y^/yRT 


(5) 


where y is the ratio of specific heats. This choice to obtain 
simplification of the equations is not unique; similar advantages 
result with other dimensionless combinations of velocity squared and 
a temperature. For instance, some workers have \ised the ratio of 
dynamic temperature to total 'cemporature; in appendix B of the 
present report are presented the canonical differential equations 
in terms of this variable. 
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If equations (1), (2), (3), and (4) are divided by CpT, pVA, 
p, and p/EpT, respectively, there result — 


xL. — ^ 

yRT V T " CpT 


dY dp dA 

dV dp pdF 

^ yET V ^ p p 


dp dT dp 

p T p 


( 6 ) 


( 7 ) 

( 8 ) 
( 9 ) 


With use of equation (5) and the expression for dV/V obtained by 
logarlthmio differentiation of eqiAatlon (5), 


dV l^dJN dT\ 

V 2 V N T / 

and, upon elimination of dp/p, there are found 


( 10 ) 


(7-1)N dN 


1 +. 


(r-i)N 


dT dQ 


= dS 


( 11 ) 


N L 
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J 

1 T 

CpT 


1 dN 

dp 


L dT 
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— 

+ — 

- - 
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- — =; da 

(12) 
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7N dN 

dp 

vS dT 

dF 


« - 




— s dp 

(13) 

2 N 

P 


T 

ET 



where the dimensionless quantities d9, da, and dp, have been 
introduced to simplify the following analysis. 


n 


h 
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Solution for Logarithmic Differentials 

If the determinant formed hy the coefficients of AN/n, dp/p, 
and dl/T in equations (11), (12), and (13) is not identically zero, 
the equations may h© solved uniquely for these, three differentials. 

As the determinant in question is proportional to (l-N), -which vanishes 
only for N = 1, the solution is obtained as follows: 

dN/N = ( 1 -N)~^ f(l + 7 N) d 0 + [2 + (/ - 1 )n] dn + [2 +( 7 - 1 )n] da^ (14) 

^p/p = (1“N) ■<1 - 7 N d 6 - 1^1 +( 7 - 1 )n]- dp - 7 N da^ (15) 

dT/T = (1-N)"^ I (I- 7 N) dS - ( 7 - 1 )N dp - ( 7 - 1 )N da ^ (16) 

It is also convenient to record the differential expressions for 
the density p and velocity T: 

dp/p = dp/p - d£C/D = (1-N)“^ (-de - dp - N da) (1?) 

dV/V = (dN/N + ar/T)/2 = (l-N)"^ (d0 + dp + da) (18) 

Application of Second Law of Thermodynamics 

The first law of thermodynamics -was used in the formulation of 
the basic equations; the second law of thermodynamics may be 
employed to furnish additional information. The entropy differ- 
ential dS for a perfect gas is given (cf. reference 13 , p. 63) by 

dS/Cp = dT/T - ( 7 - 1 ) /r] ^P/P = d0 + [(7-l)/7j (19) 

The secpnd law of thermodynamics then states . . 

0^ dS/Cp - dQ/cpT = [(7-1) /r] ap (20) 

The relation, according to equation (19), that 

dS/Cp = (d 0 + dp + da) - dp/y -da 
when used with equation ( 20 ), results in the inequalities 
d 0 = dS/c^ = dp - dp /7 - da = dp -da 
where dp = - d 0 + dp + da. 



( 21 ) 
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DISCUSSION OF EQUATIONS 
Remai’ks on Integration of EcLuations 
Equations (14), (15), and (16) can be rewritten as 


N' 

N 

l+yN 
— ^ Q' 

N 

2+(r-i)N 

N 2+(7-1)N 

(22) 



CpT 


ET 

1-N A 

P' 

1 

L2?L Q, 

_P_ 

1+(7-1)N ^ 

• A« 

(23) 


1- 

■N c T 
p 

1-N 

ET 

1-N A 

T' 

T 

= T::w 

1-yN 

Q' 

c^T ^ 

T 

1-N 

(r-i)N 
ET ^ - 

T (r-l)N 
1-N A * 

(24) 


where the primes indicate differentiation with respect to x. This 
system clearly aatiafles, except at N = 1, the conditions of the 
fundamental existence theorem (see, for example, reference 14) when 
Q, F, and A are differentiable. Hence a solution exists except 
at sonic velocity and may be obtained formally wheii possible, and 
by standard n;mnerical methods otherwise, as soon as the functions 
Q, F, and A (or their derlvati-res) arc specified. More gener- 
ally, the system may bo solved in similar fashion for any thi’oe of 
the variables N, p, T, Q, P, and A as functions of x, when 
the variation with x of the other three is prescribed. Also it 
may be noted that as all the foregoing variables are functions of 
one parameter, etny two may be considered as functions of each other 
under siil table circumstances. 


Direction of Change of Flow Yarlables 

In practical as well as in theoretical work it is frequently 
useful to be able to determine the direction of change of flow 
quantities with respect to heat addition, friction, or area varia- 
tion without troubling to get quantitative information from inte- 
grated forma. Equations (14) to (18) (or (22) to (24)) permit the 
specification of signs of derivatives at any particular point and 
also throughout certain regions of flow. Thus equation (14) shows 
that in subsonic flow the effect of positive 0', p', or a', 

is to increase N, whereas for supersonic flow the effect is to 
decrease N. When the derivatives have different signs, the net 
effect will depend upon the algebraic sum of the separate con- 
tributions. 


I 
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As an example of tlie use of this technique, suppose heat is 
added to a fluid in a consteint -area pipe, with negllgihle friction; 
that is, 6 ' 3^ 0, |i’ =« a' = 0. It is easily seen from equa- 

tions (14) to (18) that for the entire range of N from zero 
to infinity 


(1-N) dN/dQ ^ 0 
dp/dN = 0 
(1-vN)"^ dT/clN = 0 
dp/dN = 0 
dV/dlf = 0 


(25) 

(26) 

(27) 

(28) 
- (29) 


These results are given in reference 12. By use of the chain rule 
for the derivative of a function of a function, the sign of the 
derivative of any of the flow vai'iahles with respect to any of the 
others may he obtained; thus, from equations (25) and (29) it is 
clear that 


. , d7 . , dY dN 

(i-“) Sq * 


(30) 


As another example, consider the flow in circular cylindrical 
pipes with heat addition and with friction; that is, 6 ’ ^ 0, 
li' 3^ 0, a' = 0, (See also related discussion in reference 8.) 

Equation (14) will he used to determine the direction of change of 
N with respect to x. If the heat addition is only through the 
wall, which is at temperature T^, the heat added per unit mass of 
fluid in passing a distance dx along the tube is given by 

PYA dQ = h (T^-T) (rtD dx) (31) 

where D is the tube diameter, and h the local surface -to -fluid 
coefficient of heat transfer, in heat unite transferred per unit 
temperature difference, per unit area. (Cf. equation (2), refer- 
ence 15, p. 135.) In conjunction with equation (11), equation (31) 
leads to 


- 7(D dx _ 4h C(T^/t) - I] dx 

CpT Cp pYitD^/4 Cp pYD 


( 32 ) 
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The expression for frictional work done is assumed to be elven (cf . 
reference 15, p. 119, equation (8)) by 


fV^ dx 

dF = 

2 CD/4) 


(33) 


where f is the Panning friction factor. From equations (13) and 
(33) it follows that 


dF 2fYN dx (-zas 

^ ^ 

If Reynolds' analogy la valid, h may be replaced by Cp pTf/2 
(reference 15, p. 162, equation (1)), whence equation (14) becomes 

dN/dx =Cn/(1-N)] |^(1+7N) [ (^v/t) - l] + [2»- (r-l)N j yN ^2f/0 (35) 

This equation may be used to determine the direction of change of 
N with X, and hence of other flow quantities, for various 
ranges of N and of T^/t. For values of (T-^f/T) «1 (maximum rate 
of cooling), dN/dx is positive for values of 

1> N > j^-7 j/27(r-l) =.0.58 

for y = 1.4; that is, the effects of friction in increasing the 
Mach number over'^lance the effects of the colcU walls in lowering 
it if 1 >M = vN> 0.76 for y = 1.4. If N > 1 then (dN/dx) < 0, 
and acceleration of frictional, supersonic flow by convective 
cooling appears to be impossible. Acceleration of frictionloss 
supersonic flow by cooling should, however, be possible (refer- 
ence 12 ) . 


Behavior of Solution at Sonic Velocity 

The differential equations (14) to (18) must be examined 
as to behavior at the singular point N = 1. In order that the 
logarithmic differentials may be defined at this point, it is neces- 
sary that dp = de + dp + da veinish suitably at N = 1; that is, 

dp =de + dp + dix = 0 at N = 1 (36) 

because each logarithmic differential is proportional to dp/(l-N) 
there. If dp 9 ^ 0 upstream of the end of the duct, N can equal 


T 
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1 only at the end of the duct. This situation le illustrated hy 
the "choked" converging nozzle and by the frictional diabatic floWj 
which is treated in the previous section. Eq.uation (36) is formally 
satisfied at the end of a duct where d6, d[i, and da may be 
considered to vanish for all values of N. 

Between the ends of a duct^ however, dP must always vanish 
where N = 1. This condition shows that at N = 1 arbitrary varia- 
tions of d©, d|j., and da are not possible. A specific illustra- 

tion is the ideal nozzle in which d© = d^i = 0; according to equa- 
tion (36), da is then restricted to the value 0, which means 
that the area has a stationary value at N = 1. This is the well- 
known result that sonic velocity can be attained only in the throat 
of an ideal nozzle in shockless flow. A quite similar treatment 
applies for the cases where d© and dja are the quantities to be 
investigated . (See pertinent material in references 5 to 12 . ) 
Condition (36) , which was necessitated by the presence of the deter- 
minant (1-N)/2 in equations (14) to (16) is thus seen to 
provide a unification of the treatment of the flow behavior in the 
neighborhood of sonic velocity. 

Combination of the second law of thermodynamics with equa- 
tion (36) also yields limitations on the behavior of the flow at 
N = 1. According to equations (21) and (36), at sonic velocity 

d© ^ dS/cp § - da = dA/A (37) 

These results may be stated in words to the effect that in converg- 
ing or constant-area channels at N = 1, neither the heat term 
d© = dQ/cpT nor the entropy term dS/cp can be positive. In 
I diverging channels these two terms may be either positive or neg- 
latlve. If either d© or dp, is everywhere 0, relation (21) 
yields more detailed results. For example, if d© s 0, then at 
N = 1, by equations (36) and (20) 

dp = -da 

= [7/(7-!)] dS/cp ^0 (d© = 0, N = 1) 

Continuous flow with friction and without heat addition at sonic 
velocity cannot therefore occur in a converging channel (refer- 
ence 7) . 
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A more complete treatment of the hehavior of the flow when 
N = 1. between the ends of the duct is obtained by considering 
second-order terms. As N approaches unity, eq.uation (14), which 
may be written 

W’/N = j(l+7N)0' + [^+(7-1)N] u' + [2+(7-1)N] a' l/Cl-N) 

(where the prime indicates differentiation with respect to ^), 
takes the form o/o. For the evaluation of this limit, L’Hopital's 
rule gives after some calculation 

Wo’ = [©o' No' + (1+7) Po"3 /(-No') 

where subscript o denotes value of function at N = 1. The solu- 
tion for Nq ' is 

No ' = -0o ' /2 ± y<0o72)^ - (1+7)3 o" (38) 

The double-valuedness of the derivative at N = 1 wilL have impor- 
tant consequences in that a unique solution of the equations may 
not be obtained when N » 1 along the flow path. In general, it 
will be possible to continue the solution from N = 1 along either 
of two paths, depending on the choice of sign. In certain 'cases, 
depending on the signs of $o' Po"^ sign will correspond 

to continuation into subsonic flow, the other into supersonic; 
otherwise the two choices will correspond to different continuations 
into flow of the same character. This result means that specifica- 
tion of initial conditions and of variation of d0, dp, and da 
alone is not suffio§g,nt to insure a unique solution if N becomes 
unity along the flow. In the event that the radicand is zero, it is 
possible that only one solution is obtained; or it may happen that 
some higher derivative is double-valued with resulting ambiguity of— 
solution. The analysis for this case is somewhat Involved and will 
not be continued here. 

It is interesting to note that a less general problem of the 
same nature has been presented by Lorenz (reference 16) and Prandtl 
and Proell (reference 17) . Some of this work is possibly more 
accessible in reference 18. 


The Phenomenon of Choking 

The general equations (ll), (12), and (13) Impose restrictions 
on the relations between the flow variables and the heat, friction, 
and area variation. When these restrictions take the form of upper 
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or lower limits on mss flow, the associated phenomena are termed 
"choking" processes. As an example, it is well known that the ideal 
nozzle has for given suhsonic entry conditions a maximum mass flow 
beyond which the discharge cannot be increased no matter how much 
the exit pressure is lowered. Another case is "thermal choking, " 
wherein the entrance Mach number and hence mass flow in diabatic, 
frictionless, constant-area flow is limited for given heat addition 
despite indefinite reductions in outlet pressure (reference 12). 

The nature of choking may be studied with the help of eq.ua- 
tion (22), which was derived simply from the basic eq.uations. It 
will be shown that unless heat, friction, and area variation are 
such that (1-N) times the right-hand side of eq.uation (22) 
changes from positive to negative as x increases, the Mach number 
in the tube cannot become greater than 1 if the entrance velocity 
is subsonic and cannot become lees than 1 if the entrance velocity 
is supersonic, provided that the flow variables remain continuous. 

For convenience, designate by Y the factor (l-N) times the 
right-hand side of eq.uation (22). The q.uantity Y is seen to con- 
sist of a sum of terms in Q', F', and -A' multiplied by func- 

tions of N that are always positive. (In the event that only one 
of the te 3 mis Q', F', and -A' is not 0, Y becomes merely the 

derivative of the heat added, the frictional work, or the area, 
multiplied by a simple function of the flow variables; then posi- 
tive Y corresponds to the case of heat addition, friction, or a 
converging duct . ) 

Suppose first that Y is always negative. Then if the flow 
at the entrance section x^ is subsonic, dN/dx = (1-N)Y<0, and 
the Mach number decreases; if the entering flow is supersonic, 
dN/dx = (1-N)Y >0, and the Mach number increases. 

Suppose now that Y is always positive. Then, if the entering 
flow at X]_ is subsonic, dN/dx = (1-N)Y>0, and the Mach number 
increases. But N cannot increase past unity as x increases. 

For suppose N = 1 at x = Xg and is greater than 1 in the right- 
hand neighborhood of Xg (exclusive of Xg); then dN/dx is nega- 
tive in this neighborhood, because (l-N) is less than 0 and Y is 
greater than 0. Now N is eq.uai to unity at x = Xg, is contin- 
uous, and has a negative derivative in the neighborhood mentioned. 
Hence N is less than 1 in this neighborhood, which contradicts 
the assumption. Therefore N cannot be greater than 1 if Y is 
always positive, N is continuous, and N(xj^) is less than 1. In 
general, no continuous solution exists for values of x > Xg if Y 
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ia always positive. This statesnerit , and the foregoing proof, aro 
valid ovon if N'(Xq) does not exist. Au analogous development 
may bo made for W(X]^) greater than 1 with the conclusion that, 
with Y positive and N continuous, N cannot bo less than 1. 

If Y changes from positive to negative at x = Xq, however, 
the value of N may cross unity at that point, but if Y is ini- 
tially negative, K goes away from unity as previously shown and 
can only turn toward xmity if Y changes from negative to positive; 
this change must he made at some value of N other than 1. After 

Y has changed to positive, the situation reduces to tho case that 

Y is always positive, with tho N where Y changes sign now taken 
as the entrance N. 

It has been shown that up to some fixed point in tho tube, 
which can bo either the exit or the point at which Y changes from 
positive to negative, N and -therefore the Mach number do not 
become greater than 1 if tho entoi-ing velocity is subsonic nor loss 
than 1 if tho entering volcoity ia supersonic. Furthermore, if Y 
la positive up to the point at which N is limited, the derivative 
of N before this point is always positive if tho entering flow ia 
subsonic eind is always negative if the ontering flow is supersonic. 
Thus, for posltivo Y and subsonic entrance velocity the entrance 
N cannot exceed some llmitirg value less than 1 determined hy tho 
particular Q, P, -A variation, for N is alwajs increasiiJg 
from its initial jsaluo and cannot exceed 1; eind, by analogous con- 
siderations for positive Y and supersonic entrance velocity, N 
cannot be less than some limiting value greater than 1. This 
limitation is essentrlally the choking phenomenon. 

The specific form this limitation takes is not easily stated 
in the general case, because the choice of wiilch factors aro to.be 
held constant and of which vnrlahlos are to be considered limited 
determines the particular fonn of the rostnictions. Numerical 
results for some -particular cases are given in reference 12, among 
others, to illustrate tho nature of possible results. It is felt 
that additional special cases should be investigatci before a thorough 
study of the general case is attempted. 


Flight Propulsion Bosearcb Laboratory, . 

National Advisory Committ^ for Aeronautics 
Cleveland, Ohio, February 24, 194?. 
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APrENDIX A 


DEEXVATIOWS OF THSEE BASIC EQUATI0R3 
Conservaticfn of Energy 

The first lav of thermodynamiCB applies to energy changes 
between two states of a systsEi enclosed within a surface* Let the 
systeia (fig. 1) be the gas of mass Am. that is contained in the 


X = 0 



Initial state within the tube walls and the sectione at end 

^2 and in the final state within the tube walls and the sections 
at and ^1^ ^1 arbitrary, and 

determined by the condition that the mass between and t 

equal the mass between Xg t 2' ™(^) is defined as the 

total mass of gas contained within the duct between the sections 
X = 0 and x = x, the definition of Am becomes 

Am=m(x 2 ) - m(x^) ^ :3lH (33) 

Let U(x) denote the internal (thermal) energy' per unit mass 
of fluid, p(x) the static pressure, V(x) the gas velocity. 
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A(x) tb.9 croes-sectional area of tha duct, each.-tai:en at the .sec- 
tion xj and Q(z) the heat ( in mechanioal-enorgy units) added 
from walls or by combustion^ to a unit mass of the fluid during its 
passage from x = 0 to x = x. Then the first law of thermodynamics 
says for steady flow that 



As A dx = (l/p) dm, where p(x) is the density of the gas, equa- 
tion (40) becdmes 



(41) 


when the limits of integration are changed. 

Provided that the integrand is continuous (which requirement 
excludes shock), the expression on the right-haiid side of equa- 
tion (41) may be written, by the theorem of the mean for integrals: 

- m(x3^)J f(ni_L*) 

where f (m^^*) indicates the value of the integrand at some m = Dii*^ 

^(^i) ■*^®1*'^ l) . The integral on the left-hand side yields a 

similar result, with subscript 1 replaced by subscript 2, whence, 
by virtue of equation (39) 

f(mi*) = f(m2*); m(xi) < mi* m(^l), m(x2) < O-z* <■ “(^g) (^2) 




Actually the heat liberated by cambustion might be considered as 
part of internal energy; or the external surroundings ml^t be 
considered to Include the fuel; or the first law mlgfit be genera- 
lized to include heat sources. The treatment given hero is con- 
venient but must be understood to require some Justification on 
one of the bases mentioned. 


r 


I 


k 
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As I 2 — ^^2^ this eq.uation beccmes 

f = f^Cxg)] (43) 

but z-j_ and. ^2 ®-xe arbitrary, HenCe f is a constant, and 

j|(u + 1 v2 . e + p/^ . 0 

or, since d(U + p/p) = dE, ■Hhere H is the enthalpy per unit 
mass 

dQ = dJE + VdV (44) 

For a perfect gas, dE = CpdT, "whence the energy equation is finally 

dQ = o^dH + VdV (45) 

Conservation of Mass 

The conservation of mass, in the form useful here, states that 
in steady flow the mass entering a closed surface during any time 
Interval At is equal to the mass leaving during the interval At. 
Let the closed surface consist of the sectlans at arbitrary | ^ and 
X 2 (fig. 1) and the portions of the duct, between these sections and 
let At be the time required for the mass m^^^) - m(i^) to enter 
the surface while mass "^^2^ ~ flows out. The value of 

is arbitrary, whereas ^2 is fixed by the conditions on the time 

Interval. Upon definition of t(z) as the time required for a 
fluid particle to travel from origin x = 0 to z = i. At may be 
defined by 

At = t(|j^) - t(z^) = t(|g) - t(Zg) (4:6) 

The law of conservation of mass says that 



Upon change of variable, this equation beccm.es 
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A ^ /14. 

PA dt 


(48) 


as before, by the theorem, of the mean for integrals and. condi- 
tion (46) the integrand must be constant, whence 

d(pVA) =0 (49) 


Equation- of Motion 


The vector form of the second lav of motion for continuous 
media states that the integral of the density of surface forces 
over a closed s\irface is equal to the integreil of the density times 
the particle derlvativs of the velocity over the volume enclosed by 
the surface, (See, for example, reference IS.) For the mass of 
gas contained vithin the sections at x^_ eind Xg (^i^* 1) 
walls of the duct, the horizontal component of the equation of 
motion becomes 


p(xi) A(xj_) 



(50) 


where B(x) dx is the force exerted on the gas by the portion of 
the duct be-tween x and x + dx. For steady flow, the integral on 
the right-hand side may be transformed as follows: 





pVA 


■r— dx 
dx 


= pVA [t(x2)-T(Xi)J 
(51) 


It follows from equations (5p) and (51) that the equation of motion 
may be written in differential form as 

pdA + Adp + pVAdT = Edx (52) 


Wow Edx may be resolved into two constituent parts, the 
horizonteil components of the tangential friction drag and of the 
normal pressure reaction. If the half-emgle of the duct is denoted 
by 0(x), the wall surface by a(x), and the tangential frictional 
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drag per unit area ty T (i), from figure 2 it is clear that 



and that 


da = dA/sin 0 


hence 


R dx = -(t do) COB 0 + p dA (53) 

It is poBsihle to use equation (53) directly without further 
analysis if the friction factor is related to T hy the equation 

T = fpV^/2. In many engineering treatments, however, f is defined 
in tems of the "energy loss due to friction." In order to make 
this concept of energy loss more precise and to make possible 
derivation of a rigorous connection between t and energy lose 
define F(x) as the work done by unit mass of the fluid against 
friction in moving from the origin to position x. The work done 
in moving the entire mass of fluid between Xj_ 
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and ^ to tbe region betveen X 2 and.-Js (fiS* 3) vill "be cesa- 
puted in torme of the original variables T and ex' and in terns 



Figiiro 3 


of the new variables F and m. If the two quantities are equated 
and suitably transf orned ^ a relation will be obtained between dF 
and TdO. 

Let ^ 2 _) and ig be chosen arbitrarily^ and let be 

determined by the condition that the toteil mass Am between 

1 1 equal the total mass between Xg and $ g . Let Xg^ be an 
arbitrary point between Xj_ and 1 3 _, and let be dotemlned 

by the condition that m(Xg) - m(xjL) = m(x^) - m(xg). In particular^ 
if_.x^ = then x-j^ = Xg; if Xg^ = then x-^j =£g. Thus it 

is sden that^ as x^ runs from to x-j^ runs fraa. ig to 

^ 2 - 

In order to determine the Trork done in terms of the variables 
T and a., the procedure is to move thin sections fresn their original 
positions, given in each case by to their final positions 

to find the work done by each of them, and then to add up the work 
done for all the sections. 
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First ^ divide the mass Am. hstveen Z 3 _ and into smaller 

elemients of mass • • • ^n®'* (^®® fig- 3.) Let A^a 

be the wall surface corresponding to the element of mass A^m and 
let the duct wall between the initial and the final position of 
Aj^m be divided into elements of length A<s. Fop a given A^m it 
follows from figure 3 that 

A^m = Aj_ 

A.<y = 2jtr. A.s 

X 11 

A^i = cos A^s 


so that 

■A^0 = g^ A^m 

where giS. Zitv^/p^ cos 0^. Each of the quantities 
and gj_ is to be evaluated at the proper point within Aj^x. 

Finally from figure 3 there exists on A^n a force Ao'j_(sj), 
to be overccme by the work against friction, equal to 


t(Sj) A^(sj) =t(bj) gi(Sj) A^m 


where it is convenient to consider t and o as functions of s, 
inasmuch as frictional dreig and the extension of an element of fixed 
mass depends on the position of the element in the tubes. 


For anj' element of mass Aj|^m then, an approximation (the 
acciuracy of which depends on the size of Aj^m) to the work done 
by A^m in going frcxa section Xg^ to section x-jj is 


11m 
k - > oo 
A^s-^0 


s(Xb) 

/ = 

s(Xa) 


lim k 

k >09 

A^s->0 



t(Sj) g^Csj) Aj^m AjS 


= Aj^m 



t(b) g^Cs) 



da 
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It follows that the oorresporidirig, approxlmaticii to the work dono 
by the entire mass Am is 


n 


i=a 


2 /t(o) gj(B) de>4jm 


and, as- the approximation becomes mope and more accurate as the 
largest Aj^m becomea smaller and smaller, the expression for the 
work, done against friction when tlve entire mass between and 

is transported to the position between 13 and ig becames 


lim 

n— 

Aj^m — >0 


^(^ 1 ) f s[®(^b)J 


/t(8) e,(a) da > Aim =| I 7 ( 0 ) g(0) ds dm 


M^i) 


The limit'“of the double sum may be' written as an iterated integral 
as follows : 


mUi) / s[m(x^ 5 )] 
U(xi) Js[m(xa)] 


Tg ds dm 


(54) 


It is clear that the work done against friction when the entire 
mass is moved as previously described is equal also to 


11m ^ ^ 

AjF— 3>0 y y AjF Aj^m = 

A^^m-^O 1 J 


r F[m(rij)] 




dF dm = 


F[m(x|jJ] 


'm(xt) -^F[m(x )] 


'm(xi) JF[m(xa)] 
(55) 


ds 


dsdi 


whence it follows, since the Intervals for both integrations are 
arbitrary, that 
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and 


cLF* 


/ 2Tcr \ / dg \ 
*^\pA COB j \2nrj 


then 


(xda) COB = coa^ pA dF (57) 

Because, fpr the small angles usually under consideration, cos^ ^ 
is very nearly unity (for a half angle of 6°, cos^ ^ = 0,989), the 
retention of this factor except for particularly precise work would 
not seem Justifiahle. Hence equation (53) may he written 


E dx = -pA dF + p dA (58) 

The differential equation of motion is finally 

-dp = pVdT + pdF . . (59) 

The connection between dF and the differential loss in 
stagnation pressure (-dpt) can, with the help of equations (14) 
and (15), be expressed in the form 


^ + -rm (so) 

Pt RT [2 + (7-1 )n3 

Thus except for the limiting case of incompressible flowj (-dp^) 
and dF cannot be used interchangeably in defining the friction 
factor even for the adiabatic case ■»diere dS = 0. 
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APPENDIX B 
THE Z LANGUA13E 

In the place of - N = / 7 ET the equations may be formulated 
in terms of Z, defined through 


Z = (v2/2cp)/(T + ’V^/2Cp) 


(61) 


The numerator is the so-called dynamic temperature, the 
is the total temperature. Z and N are related by 

denominator 


z = [(r-i)N]/C2-^(r-i)N] 

(62) 


N = 2Z/(7-1) (1-Z) 

(63) 



This replacement represents a one-to-one transformation of N into 
Z in the range 0 to infinity for N, 0 to 1 for Z. 

In order to illustrate the form that some of the earlier equa- 
tions take under this transformation, equations (14), (15), aixd (IG) 
are written in terms of Z, 



(64) 


(65) 


dT 


2 



1-Z tZN 

7 d© - Z du 

2 7-l-^ 


- Z da 


( 66 ) 


T 
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